Finite-Size Scaling Analysis of the Eigenstate Thermalization Hypothesis in a 

One-Dimensional Interacting Bose gas 
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By calculating correlation functions of the Lieb-Liniger model based on the algebraic Bethe ansatz 
method, we conduct a finite-size scaling analysis of the eigenstate thermalization hypothesis (ETH) 
which is considered to be a possible mechanism of thermalization for isolated quantum systems, and 
find the ETH holds in the thermodynamic limit even for an integrable system. However, we show 
that the ETH gives only a logarithmic correction to other thermalization mechanisms such as the 
eigenstate randomization hypothesis and the typicality argument for the microcanonical ensemble. 
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Introduction. — Recently, there has been a resurgence 
of interest in understanding thermahzation from quan- 
tum mechanics due in part to highly developed experi- 
mental techniques in ultra-cold atomic gases which en- 
able one to prepare, control, and measure isolated quan- 
tum systems [l|, B, and also due to several theoretical 
breakthroughs [3|-[l0[. Among the underlying mecha- 
nisms proposed so far for thermalization in isolated quan- 
tum systems are the eigenstate thermalization hypothe- 
sis (ETH) [ll-13L the typicality of quantum states with 
large dimensions [y, MK and the eigenstate randomization 
hypothesis (ERH) [m^^ The ETH, which has attracted 
a growing attention |lll4l3j , states that the expectation 
value of an observable stays the same over many-body 
eigenstates having close eigenenergies. The validity of the 
ETH has been examined by using the exact diagonaliza- 
tion of the Hamiltonian 13 - 13| • It has been claimed that 
the ETH holds for the case in which the system is non- 
integrable or chaotic 13|, llSMlTf . These studies focus on 
the dependence on the characteristics of the system (such 
as geometrical configurations, sets of parameters, etc.) 
for a fixed size of the system. On the other hand, the 
ERH states that the expectation value of an observable 
over consecutive many-body eigenstates randomly fluc- 
tuates, and the typicality argument concerns statistical- 
mechanical properties that are shared by randomly gen- 
erated quantum states from a large-dimensional Hibert 
space. Although all of the ETH, ERH, and typicality 
argument can explain thermalization, it has not been 
quantitatively discussed which of those mechanisms dom- 
inantly contributes to thermalization. To resolve this 
problem, it is essential to analyze the finite-size scaling 
of the ETH. However, this goal has not been achieved 
yet because the numerical cost to diagonalize the Hamil- 
tonian grows exponentially with the size of the system. 

In this Letter, we identify the finite-size scaling prop- 
erty of the ETH by calculating the correlation functions 
at finite temperatures for the Lieb-Liniger model 18 1 



which allows us to obtain the exact many-body eigen- 
states by the Bethe ansatz and thus to circumvent the nu- 
merical difficulty in diagonalizing the Hamiltonian. We 



show that the ETH becomes better satisfied for larger 
systems. This suggests that the microcanonical ensemble 
is still useful to calculate non-conserving quantities in the 
thermodynamic limit, while the generalized microcanon- 
ical or Gibbs ensemble lead to better predictions in small 
integrable systems 3"21|. Then, by using a finite-size 



scaling of the ETH, we discuss the quantitative relation 
between the ETH and other scenarios of thermalization 
in isolated quantum systems such as the ERH and 
typicality argument [1, Q • We find that the ETH con- 
tributes to thermalization at most as a logarithmic cor- 
rection to the ERH and typicality. 

Model. — The Lieb-Liniger model describes a one- 
dimensional Bose gas with a delta-function interaction. 
The Hamiltonian is given in units oi h = 2m = 1 as 



H= i dx[d^-^\x)d^-^{x)+c^\x)-^\x)-^{x)-^{x)\, 
Jo 

(1) 

where ^'(a::) is the bosonic field operator, L the linear di- 
mension of the system, c the magnitude of the contact 
interaction, and the periodic boundary condition is im- 
posed. The iV-body eigenstates are constructed from the 
monodromy matrix [23] given by 



A{X) BiX) 
C{X) D{X) 

= : V exp 



da; 



A/2 ^c^^{x] 
-^/^^{x) -A/2 
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where V and : • • • : denote path ordering and normal 
ordering of bosonic field operators, respectively. The 
TV-body eigenstate can be obtained by acting B{X) N- 
times on the Fock vacuum |0), |{Aj}) = Wj^i B{Xj) |0), 
provided that the set of parameters {Xj}jLi satisfy the 
Bethe equations 
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(j = l,2,...,iV). (3) 
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Since the energy and momentuni of this state are given 
as E = J2f=i^^ and P = J2f=i^j^ respectively, the set 



of parameters {Xj}jLi may be regarded as the momenta 
of the dressed noninteracting particles. 

Correlation Function. — We consider the expecta- 
tion values of the quantity ^''''(x)5'(0) over each eigen- 
state, i.e., ({Aj}|^'1'(x)*(0)|{Aj}). Substituting 'i'^x) = 
e~'^^5'(0)e'^^ and inserting the completeness relation for 
the {N — l)-body eigenstates, we obtain 

({A,}|vI/t(^)vi,(0)|{A,}) 

{Mfo }jv-l 



(4) 



where the summation is taken over all the {N — 1)- 
body eigenstates. The form factor ({Aifc}|^'(0)|{Aj}) 
is reduced to N scalar products {{f^k}\{^j}j^i) (l = 
1,2, . . . , N) which are represented by the determi- 
nants of {N — 1) X {N — 1) matrices [25[- The crucial 
observation is that these determinants can be summed 
up, giving a single determinant of an x matrix (26| . 
This algebraic manipulation greatly reduces the compu- 
tational task, enabling us to conduct a finite-size scaling 
up to a decent value of N. 

By increasing N and L with their ratio held fixed, we 
calculate the real and imaginary parts of the quantities 
({Aj}|^'t(a;)*(0)|{Aj}), as plotted in Fig. [H We set the 
parameters as c = 10 and x = 1/2, where the unit of 
length is taken to be the mean distance between parti- 
cles L/N. The energy window is taken as [Eg, Eg + 10], 
where Eg is the ground-state energy which depends on 
the size of the system. We have plotted all energy eigen- 
states in the energy window. Figure [T] clearly shows that 
the fluctuations of the eigenstate expectation values be- 
come smaller with increasing the size of the system. Since 
the ETH states that such fluctuations vanish in the ther- 
modynamic limit, we will refer to these fluctuations as 
the ETH noise, which is analyzed as follows. First, since 
the real part exhibits a systematic linear behavior (see 
Fig. [T^), we first identify the linear part by the least 
squares fitting and subtract it from the data. Then, the 
root mean square of the subtracted data gives a quanti- 
tative measure of the ETH noise which we will denote as 
a. The dependence of the ETH noise on the size of the 
system is illustrated in Fig. [5^. This figure shows that 
the ETH noise decays with the number of the particles 
as (T oc N~", where the exponent a is obtained through 
best fitting as a = 2.06 ± 0.01. Second, we evaluate the 
ETH noise for the imaginary part from the root mean 
square of the data (see Fig. Wp). As with the real part, 
the ETH noise decays with the number of the particles 
as a power law (see Fig. [Hd) with a = 0.784 ± 0.006. 

These results for the real and imaginary parts of the 
eigenstate expectation values of the quantity ^'^{x)'i'{0) 
shows that the ETH becomes better satisfied as we ap- 
proach the thermodynamic limit despite the fact that the 




o 

1? 



0.4 
0.2 


-0.2 
-0.4 
-0.6 




N=5 
N=10 X 

N=15 + 
_i I I I 



10 15 20 25 30 35 40 

Eigenenergy [?i^7V^/(2mL^)] 



45 



FIG. 1: (Color online) (a) Real and (b) imaginary parts of the 
expectation values of ^'^(a;)^(0) plotted against the eigenen- 
ergy for N = 5 (asterisk), 10 (cross), and 15 (plus), where 
X — L/2N. For each case, all the energy eigenstates in the en- 
ergy window [Eg, Eg -\- 10] are shown, where _Eg is the ground- 
state energy. The fluctuations of the data around the fitting 
lines, which we call the ETH noise, become smaller as the size 
of the system becomes larger. 
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FIG. 2: (Color online) ETH noise plotted for real (plus) and 
imaginary (cross) parts of the correlation function as func- 
tions of the size of the system between N = 5 and — 15. 
For both the real and imaginary parts, the ETH noise de- 
cays as (J oc N~°'. The exponents are obtained by fitting as 
a = 2.06 ± 0.01 and 0.784 ± 0.006 for the real and imaginary 
parts, respectively. 
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Lieb-Liniger model is integrable, whereas previous stud- 
ies [l3 , 15 - 17 1 showed that the ETH holds worse in inte- 
grable systems than in non-integrable chaotic systems at 
a fixed size of the system. 

Interplay between ETH and ERH. — We now discuss 
which of the ETH and the eigenstate randomization hy- 
pothesis (ERH) [22'] dominantly contributes to thermal- 
ization. The ERH states that the expectation values 
of observables over each energy eigenstate fluctuate ran- 
domly. If the ERH holds, the following inequality holds 
which gives an upper bound on the difference between the 
long-time (LT) average and the microcanonical ensemble 
(ME) average of an observable A 22| : 



{A)^' {A) 



ME I 



(5) 



where a quantifies the ETH noise oiA as described above, 
D' is the smoothness of the weights of the initial state on 
each energy eigenstate, and A'state is the number of the 
eigenstates in the given energy window. The parameter 
7, which characterizes the randomness of the sequence 
{{Ei\A\E,)}fj^f'- (where \E,Ys (i = 1, . . . , A^state) are the 
energy eigenstates in the energy window), is defined as 
follows (see Ref. [l^ for details). By grouping the el- 
ements of the sequence {{Ei\A\Ei)}^^''{^'' m by m and 
representing each group by its average, we can define a 
coarse-grained sequence whose variance depends on m as 
a{mY . If a{m) decays as a{m) = aAm~'^ , we call the 
original sequence {{Ei\A\Ei)}^^*^^'' random with 7 char- 
acterizing the randomness. As shown in Fig. [3l the ERH 
holds for the Lieb-Liniger model with 7 — 0.55 ± 0.06 
and 1.3 ± 0.1 for the real and imaginary parts, respec- 
tively (see Fig. |3]) ^2?']. The right-hand side of ([5]) gives 
an upper bound on the difference between the long-time 
average and the microcanonical ensemble average and it 
is vanishingly small due to two independent factors. The 
first is the contribution that the ETH noise a a becomes 
smaller with the scaling behavior a a oc N~°' as shown 
above. The other factor N'll}^""^^ arises from the ERH 
mechanism. Now, we assume that D' — 0(1) in terms 
of the system size and note that A^statc ~ in general. 
Then, the upper bound of ([5]) has the following scaling 
behavior: 



exp 



7 



7+1 



-N -a\aN 



(6) 



where the first and second terms in the exponential orig- 
inate from the ERH and ETH, respectively. In this 
sense, the ETH contributes to thermalization at most 
as a logarithmic correction to the ERH. Although both 
the ETH and ERH can make the long-time average 
equal to the microcanocanical ensemble average, the or- 
der of magnitude is quite different. For the case of 
the scaling of the right-hand side of ([5]) be- 



^*stat< 

comes exp 



N -a\-aN 



so the ERH is still 



dominant as long as/3<7/(7 + l). 
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FIG. 3: (Color online) The variances of the coarse-grained 
sequences of the real (plus) and imaginary (cross) parts of 
the expectation values of '^^{x)'^{Q) plotted against each en- 
ergy eig enstate for = f 5 (see crosses in Fig. [T] for the real 
part [23] and Fig.[T}3 for the imaginary part). For both of the 
real and imaginary parts, the variance decays as cr(m) cx m 
where 7 = 0.55±0.06 and 1.3±0.1 for the real and imaginary 
parts, respectively. 



Rederivation from the typicality argument. — Finally, 
it is interesting to note that we can rederive the scaling 
law ([6]) from the typicality argument 0, 0, 0] ■ The typi- 
cality concerns statistical-mechanical properties that are 
shared by randomly generated quantum states. Here, we 
follow the typicality argument and show that the long- 
time average is close to the microcanonical ensemble av- 
erage for almost all initial states. We assume that the 
initial state is a linear combination of the many-body 
eigenstates in an energy window: \tp) = J2f=i*° l-^i): 
where the set of the coefficients {Ci\^^^i^° satisfy the nor- 
mahzation condition X^i^i"'" — 1- Namely, 2iVstate 
parameters Re Ci and Im Ci (i = 1, 2, . . . , A'stato) can be 
represented as a point on the (2A^state — l)-dimensional 
sphere whose radius is 1. We denote the average over 
the uniform measure on the high-dimensional sphere as 
an overline. For example, Ci = and |Cip — 1/iVstatc- 
Then, the initial-state averaging of the long-time average 
is given by 



LT 



^ \C^HE,\A\E-_ 



i=i 



^sTato E ^EM\E,) = {A) 



ME 



and the variance is obtained as 



Va = [{A) 



LT 



{A) 



ME12 



1 



(7) 



(8) 



Thus, the typical magnitude of the difference between 
the long-time average and the microcanonical ensemble 



average is given by ^/Va ~ (^aN^ 



-1/2 



Since cr^i oc " 
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and TVgtate e , the scaling of the typical magnitude of 
the difference turns out to be 

exp(^-^N - a\nN^ . (9) 

The first and second terms in the exponential correspond 
to the contributions from the typicality and ETH, respec- 
tively. This has a form similar to © and shows again 
that the ETH is a logarithmic correction to the typical- 
ity. It is remarkable that while the scaling law (|6]) has 
been obtained based on a specific model, the typicality 
argument docs not refer to the details of the system but 
only rely on the immense dimensionality of the Hilbert 
space with a uniform Haar measure. 

Conclusions and discussions. — In this Letter, we have 
conducted a finite-size scaling analysis of the ETH by 
applying the algebraic Bethe ansatz method to the Lieb- 
Liniger model, and shown that the ETH noise vanishes 
as a power law of the size of the system even for an inte- 
grable system. Thus, we have identified the quantitative 
relations between the ETH and other scenarios for ther- 
malization. We find that the contribution of the ETH to 
thermalization is at most a logarithmic correction to that 
of the ERH or typicality. This fact originates from the 
fact that the ETH is the effect of the order of some powers 
of the degree of freedom of the system, whereas the ERH 
and the typicality argument utilize the immense dimen- 
sionality of the Hilbert space which grows exponentially 
with increasing the degrees of freedom. The scaling re- 
lation ^ which includes the ETH as a small correction 
to other factors is obtained by invoking the properties 
of energy eigenstates. However, the fact that it has been 
rederived by the typicality (see Eq. ©), which relies only 
on the structure of the Hilbert space, strongly suggests 
that the obtained relation holds quite generally. This 
finding merits further study. 
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